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The surface reconstruction of SrTiO3(110) is studied with scanning tunneling microscopy and 
density functional theory (DFT) calculations. The reversible phase transition between (4x1) and 
(5x1) is controlled by adjusting the surface metal concentration [Sr] or [Ti]. Resolving the atomic 
structures of the surface, DFT calculations verify that the phase stability changes upon the chemical 
potential of Sr or Ti. Particularly, the density of oxygen vacancies is low on the thermodynamically 
stabilized SrTiO3(110) surface. 
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Perovskite oxides have attracted intensive interests in 
the fields of fundamental condensed matter physics, pho- 
tocatalysis chemistry, material science, as well as elec- 
tronics applications, due to their rich phase diagrams and 
remarkable functionalities. The related research becomes 
even more exciting since the recent discovery of a quasi- 
two-dimensional electron gas (2DEG) between two insu- 
lating materials, SrTi03 and LaAlC>3 [T]. Subsequently, 
a tremendous amount of evidence has shown that the 
perovskite oxides in the low-dimensional (LD) form such 
as interfaces, thin films, or heterostructures display an 
equally rich diversity of exotic phenomena that is related, 
but not identical to the bulk 2\ , indicating a great oppor- 
tunity for novel oxide-based devices [3] . One of the most 
important and intriguing discoveries is that the atomic 
arrangement on the surface or at the interface is determi- 
nant for the properties of the entire artificial structure. 
Hwang et al. found that the formation of 2DEG criti- 
cally depends on the type of atomic termination layer at 
the interface [T] . O vacancies (Vo 's) also sensitively influ- 
ence the density and mobility of the charge carriers at the 
heterointerface of LaA103 and SrTi03 [I]. Therefore, to 
clarify the origin of the emergent properties in LD oxides 
and ultimately to tune them for the fabrication of func- 
tionalized devices, the detailed knowledge on their micro- 
scopic structures and the high-precision growth technique 
are the key issues. 

Single crystalline SrTi03 is widely used as the epitaxial 
substrate for perovskite oxide films. In order to improve 
the growth quality, much effort has been made to obtain 
the atomically fiat and ordered SrO or Ti02 terminated 
(100) surface [SI [5]. In contrast to the electrically neu- 
tral (100) surface, the formation process of SrTiO3(110) 
surface structure is much more complicated since it is in- 
herently unstable due to the perpendicular macroscopic 
dipole formed by alternatively stacked (SrTiO) 4+ and 



(O2) 4- layers [7j. The (110) surface stoichiometry of- 
ten deviates from that in the ideal crystal, which leads to 
the formation of mixed phases of reconstruction [5] . Vo 's 
may also be responsible for the stabilization of the (110) 
surface [9 . Recently Russell et al. obtained an (nxl) 
( n=3,4,6) family of reconstructions at varying annealing 
temperatures, which was described as a homologous se- 
ries with the T1O4 tetrahedra model [TU1 [H] ■ However, it 
is still challenging to understand the stabilization mech- 
anism of SrTiO3(110) polar surface, particularly the be- 
havior of oxygen (or vacancies). 

In this letter, we investigate the microscopic structure 
of SrTiO3(110) surface and present a way of preparing 
high-quality titanate surfaces and films with ideal oxy- 
gen stoichiometry. We obtain the monophased surface 
in (4x1) or (5x1) reconstruction. The reversible phase 
transition between them is realized by tuning Ti or Sr 
concentration on the surface. With ah initio calcula- 
tions, we determine the atomic structures of both phases 
and study their stability in different chemical environ- 
ments. Thus the experimental control of surface recon- 
struction is understood within the thermodynamic pic- 
ture. Moreover, both of the experimental observations 
and theoretical calculations reveal that the Vo's density 
on SrTiO3(110) surface could be extremely low. 

The experiments were performed in an Omicron ultra- 
high vacuum (UHV) low temperature scanning tunnel- 
ing microscopy (STM) system with the base pressure of 
lxlO- 10 mbar. Nb-doped (0.7 wt%) SrTiO 3 (110) single 
crystals (12x3 mm 2 ) were purchased from Hefei KMT 
Co., China. The as-received sample was sputtered with 
Ar + beam at room temperature followed by annealing 
in UHV. Strontium and titanium were evaporated from 
Knudsen cells, respectively. The flux of Sr was cali- 
brated by observing the formation of Sr/Si(001)-(2x 1) 
reconstruction pP2] with RHEED and STM, while the 
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FIG. 1: (color online) (a) and (b) STM images (15x15 nm 2 , 
1.5 V/20 pA) of (4x1) and (5x1) reconstructed SrTiO 3 (110) 
surfaces, respectively. The insets show the RHEED patterns 
along [001] direction with (10) spots indicated by arrows, re- 
spectively, (c) The area ratio of (5x1) domains through the 
entire surface depending on the Sr dosage. The statistics are 
done over 8 STM images (500x500 nm 2 ) for each data point. 



flux of Ti was calibrated by monitoring the RHEED 
intensity oscillation during the homoepitaxial growth 
of SrTiO3(110) film [T3]. The sample was resistively 
heated by passing a direct current and the tempera- 
ture was monitored with an infrared pyrometer. DFT 
calculations were carried out with the "Vienna ab ini- 
tio simulation package "(VASP) code [TJ]. We use the 
projector augmented-wave (PAW) method [IS] and the 
Perdew-Burke-Ernzerhof functional |16j with the kinetic 
energy cutoff of 400 eV for plane waves and a (3x5x1) 
Monkhorst-Pack fc-point mesh. The surface structure 
was modeled with a supercell that was symmetrical along 
the [110] direction, consisted of a 13-layer slab separated 
by a vacuum layer of 12 A. Atoms in the centeral three 
layers were fixed and other atoms were allowed to relax 
until the force on each atom was less than 0.02 eV/A. 
Simulated STM images were generated by integrating the 
local density of states (LDOS) between Fermi level (Ep) 
and Ejr+1.5 eV in Tersoff-Hamann approximation [17] 
with the constant density method |18j . 

Figure 1 (a) shows the monophased (4x1) reconstruc- 
tion on the SrTiO3(110) surface after Ar + sputtering 
(500 eV/2 (uA for 10 minutes) followed by annealing 
at 1000°C for 1 hour. The surface appears as periodic 
stripes along the [110] direction, separated with a dark 
trench, and each contains two obvious bright rows of peri- 
odic dots. The periods along [110] and [001] are ~0.6 nm 
and 1.6 nm, respectively, corresponding to the (4x1) re- 
construction as observed clearly by RHEED. Evaporat- 
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FIG. 2: (color online) (a) The structural model of the (nxl) 
reconstruction [11] with (4x1) on the left and (5x1) on the 
right, (b) The modified structural model of (5x1). The upper 
panels are the topview while the lower panels are the sideview. 
(c), (d) and (e), (f) The simulated and experimental STM 
images of (4x1) and (5x1), respectively. 



ing a small amount of Sr metal followed by annealing at 
1000 °C, domains with a new ordering form on the sur- 
face. As shown in Fig. 1 (b), it also appears as stripes 
but wider than that in the (4x1) domain, each contain- 
ing three bright rows of periodic dots along the [110] 
direction. The periods along [110]and [001] are ~0.6 nm 
and ^1.95 nm, respectively, corresponding to (5x1) re- 
construction. The (5x1) domain area enlarges with 
the Sr evaporation dosage linearly until a monophased 
surface formed with ~0.15 ML Sr [relative to bulk- 
truncated SrTiO 3 (110), 1 ML=4.64xl0 14 atoms/cm 2 , 
see Fig. 1 (c)] . Such a phase transition can be reversed by 
evaporating ~0.15 ML Ti onto the (5x1) surface followed 
by annealing at 1000°C [T3]. 

With DFT calculations, different structural models 
of (4x1) reconstruction are compared, including the 
Sr-adatom, TiOa; and (lOO)-microfaceted models. We 
find that the TiC>4 tetrahedra model [the left part in 
Fig. 2 (a)] is energetically favorable, in consistence with 
what Enterkin et al. reported previously [II] ■ The sim- 
ulated STM images reproduce the experimental features 
very well, as seen in Fig. 2 (c) and (d). The unoccu- 
pied states STM image shows both titanium and oxygen 
atoms on the surface in bright features, as indicated by 
the ball-and-stick drawing superimposed onto the simu- 
lated image. The zig-zag TiC>4 chains along [110] corre- 
spond to the bright stripes in the STM images. 

Our experimental STM images of the (5x1) recon- 
struction cannot be described by the (nxl) homologous 



3 



<Vri < cV ) 
-4 -6 



AH„ (cV) 
-2 -4 -6 



6 


. 4x1 

5x1 




= 3 






^ 






^-3 






-(> 




(a) 



I 30c 

'J 

~. 
'- -300 
-600 



4X1 






(b) 



A^ Sr (cV) 



-8 -6 -4 
An S[ (cV) 



< " 



(C) ^ (d) 



-4 -8 -12 

A Miir (cV) 



&J*n (eV) 



-12 



FIG. 3: (color online) (a) The surface free energy of (4x1) 
(red/dashed line) and (5x1) (black/solid line) as functions of 
A/iSr and AfiTi, respectively (Afio = -2 eV). (b) The energy 
difference of (4x1) and (5x1). (c) and (d) Phase stability 
of (4x1) and (5x1) upon Sr & O and Ti & O environments, 
respectively. The black (red) zones correspond the environ- 
ment in which (4x1) is stable while the grey (yellow) zones 
correspond to that (5x1) is stable. 



series [11] very well. We modify the initial structure by 
shifting Ti04 tetrahedron 1 from the linear chain along 
[001] to the middle of the zig-zag chain along [110], as 
shown in Fig. 2 (b). After relaxation, an extra energy 
gain of 0.6 eV is obtained comparing with the total energy 
calculated for the "(nxl)"model [TJJ. Although such a 
small difference could be the result of different functional 
used in calculations, the simulated image is in excellent 
agreement with the STM observation on the monophased 
(5x1) obtained in the current work [Fig. 2 (e) and (f)]. 

In the structural models, the nominal area concen- 
tration of Ti in (4x1) is higher than that in (5x1) by 
0.1 ML. Experimentally (4x1) or (5x1) phase is formed 
by increasing the concentration of Ti or Sr, respec- 
tively. Quantitatively the evaporation dosage to induce 
the phase transition (^0.15 ML) is higher than the actual 
difference of the concentration due to different adhesion 
coefficient on surfaces of Si(100) (used for calibration) 
and SrTiO^llO). In fact the phase stability strongly de- 
pends on the metal concentration on the surface. To fur- 
ther clarify, we calculate the surface free energy of (4x 1) 
and (5x1) phases under different chemical environments, 
respectively, in terms of chemical potential. 

The surface free energy per unit area is defined by: 

7 = ^g[Gslab ~ N Ti fl Ti - N Sr [lSr ~ N Ho], (1) 

where G s i a b is the Gibbs free energy of the slab that 
can be calculated with our DFT calculations, N is the 



number of Ti, Sr, or O atoms within the slab, \i is the 
chemical potential of Ti, Sr, or O, respectively, while S is 
the area of the symmetrical surface. For the bulk SrTiC>3 
crystal, the free energy of each unit cell is: 

G h SrTi0 3 = MSr + MTi + 3/i , (2) 

which equals to the energy (Eg^ l f i03 ) that can be calcu- 
lated, since the surface is in the thermodynamic equilib- 
rium with the bulk. We use relative chemical potentials 
A/zs r , A/j,Ti and A/io, instead of the chemical poten- 
tials, with reference to the energies of a Ti atom in the 
hep bulk structure, of a Sr atom in the cubic bulk struc- 
ture and of an O atom in the O2 molecule in the gas 
phase, respectively. And the range of accessible values 
for them is determined by using the method proposed by 
Bottin et al. [15] . 

Figure 3 (a) shows the surface free energy as a func- 
tion of Afj,sr or A/iTi with A^o = —2 eV. In a Sr- 
rich (A/x Sr -> eV, or Ti-poor, A^ Tl -> -9.79 eV) 
environment, (5x1) is more stable than (4x1), while 
(4x1) is more stable in a Sr-poor (Ti-rich) environment, 
with the critical concentration Afi$r = —4.86 eV, or 
Afi Ti = -4.93 eV at A/i = -2 eV. The energy dif- 
ference of the two phases shown in Fig. 3 (b) is small. 
This explains why it has been observed experimentally 
that the two phases coexist with each other or even with 
other (nxl) phases in fragmentary domains [10 . 

At different oxygen partial pressure, the phase tran- 
sition between (4x1) and (5x1) can be induced by the 
change of Sr or Ti concentration except when A/io is ex- 
tremely low, as shown in Fig. 3 (c) and (d). At our UHV 
condition, we always observe the (5x 1) phase rather than 
(4x1) when Afj,g r is increased by evaporating Sr metal 
onto the surface. This is in agreement with the calcu- 
lation even though the phase transition is missing in an 
ideal O-free environment. However, the calculation pre- 
dicts that the (5x1) phase should be more stable than 
(4x1) when Afi^t is high and A/xq is low [the labeled 
zone in Fig. 3 (d)], totally in contradiction to the ex- 
perimental observation that increasing Ti concentration 
induces the formation of (4x1) phase in UHV. It is sug- 
gested that this is experimentally a "forbidden" zone - 
when the Ti concentration is high, the O concentration 
cannot be low. 

It is not surprising since DFT calculations have re- 
vealed that the TiC>4 tetrahedra are energetically stable 
and normally appear as building blocks of the reconstruc- 
tion on SrTiO-3 surfaces, i.e., the surface Ti tends to be 
four-coordinated by oxygen on the thermodynamically 
stable surface [TT]. We calculate the formation energy of 
Vo's on different layers, which illustrates that the energy 
on the top layer is more than 1 eV higher than on the 
buried layer. Such a characteristic is similar to that of 
V 's on anatase TiO 2 (101) surface [20] • In UHV treat- 
ment, the main oxygen source is the bulk SrTiC>3. The 
undercoordinated surface Ti atoms will bind with oxygen 
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FIG. 4: (color online) (a) An STM image (10x10 nm 2 , 2.0 V) 
of (4x1) on SrTiC>3(110). The circles indicate possible Vo's 
or hydroxyl- adsorbed Vo's. (b)-(e) The simulated STM im- 
ages of (4x1) surface with an oxygen vacancy or a hydroxy 1- 
adsorbed Vo on different sites, (f)-(h) A schematic diagram 
of the pathway of an Vo diffusing from SrTiOs(llO) surface 
towards the bulk. 



island growth when the oxygen pressure exceeds 10 -2 
mbar [23]. Moreover, specifically in MBE, the oxidation 
of the metal sources is also an important issue. The oxy- 
gen properties near SrTiC^llO ) surface region suggests 
a method to epitaxially grow high-quality titanatc films 
along [110] - as long as the metal stoichiometry is well 
controlled, Vo's can be eliminated with low growth rate 
at high temperature. 

In conclusion, we obtain atomically well defined (4x1) 
and (5x1) reconstructions on SrTiO3(110) surface. With 
ab initio calculations we resolve the microscopic struc- 
tures and explain the phase transition between the two 
phases in terms of the surface free energy determined by 
metal concentration. Moreover, the Vo's density on the 
reconstructed SrTiO3(110) surface could be extremely 
low, suggesting a feasible way to grow high-quality thin 
films of perovskite titanates. 
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atoms diffusing from bulk and form the T4O4 tetrahedra. 
Kinetically oxygen is very diffusive relative to cations in 
bulk SrTi0 3 [2T]. In Fig. 4 (a)-(c), we schematically 
show how an Vo diffuses from the (110) surface towards 
the bulk. Associated with two 3-coordinatcd Ti atoms, 
a surface Vo is unstable [Fig. 4 (a)]. It diffuses to the 
second layer and then to the fourth layer [Fig. 4 (b) and 
(c)]. The energy barrier for both processes is in the order 
of 0.6 eV as calculated by Cuong et al. [22]. Such barri- 
ers can be overcome during annealing at 1000 °C and the 
equilibrium is reached with Vo's diffusing deep into the 
bulk. 

Comparing with the theoretical simulations [Fig. 4 (e)- 
(h)], the bright spots labeled with circles in the STM 
image [Fig. 4 (d)] can be attributed to Vo's or hydroxyl- 
adsorbed Vo's (reacted with residual water) on the sur- 
face. Statistics of a number of images show that the 
possible Vo's density is < 0.4% ML, where one ML is 
defined as the number of oxygen on a perfect (4x1) sur- 
face. It is suggested that the thermodynamically stable 
surface is experimentally reached. 

It has been a great challenge to control the Vo's den- 
sity or to minimize it to obtain an oxide surface or other 
LD structures with ideal stoichiometry. Highly reactive 
oxygen sources are employed at a high partial pressure 
during the growth with molecule beam epitaxy (MBE) or 
pulsed-laser deposition (PLD) technique. However, such 
a tuning is limited since the growth mode of films changes 
from two-dimensional layer-by-layer to three-dimensional 
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